Genomic profiles of Japanese patients with vulvar squamous cell carcinoma

The incidence of vulvar carcinoma varies by race; however, it is a rare disease, and its genomic profiles remain largely unknown. This study examined the characteristics of vulvar squamous cell carcinoma (VSCC) in Japanese patients, focusing on genomic profiles and potential racial disparities. The study included two Japanese groups: the National Cancer Center Hospital (NCCH) group comprised 19 patients diagnosed between 2015 and 2023, and the Center for Cancer Genomics and Advanced Therapeutics group comprised 29 patients diagnosed between 2019 and 2022. Somatic mutations were identified by targeted or panel sequencing, and TP53 was identified as the most common mutation (52–81%), followed by HRAS (7–26%), CDKN2A (21–24%), and PIK3CA (5–10%). The mutation frequencies, except for TP53, were similar to those of Caucasian cohorts. In the NCCH group, 16 patients of HPV-independent tumors were identified by immunohistochemistry and genotyping. Univariate analysis revealed that TP53-mutated patients were associated with a poor prognosis (log-rank test, P = 0.089). Japanese VSCC mutations resembled those of Caucasian vulvar carcinomas, and TP53 mutations predicted prognosis regardless of ethnicity. The present findings suggest potential molecular-targeted therapies for select VSCC patients.

cell carcinoma in Korea 10 , which is comparable with the rates reported in Europe (69.9% in the Netherlands 11 and 69.9% in the UK 12 ).Advances in surgical techniques and chemoradiotherapy have improved the chance of survival for patients with localized disease.However, the prognosis for patients with groin recurrence, distant recurrence, or advanced VSCC has not improved in the past decades, with a reported 5-year survival rate of 25-50% 13 .Clear guidelines for the treatment of recurrent VSCC are lacking 14 , and elucidating the genomic mutational profiles of VSCC is essential to identify potential therapeutic targets.
This study investigated the genomic mutational profiles of Japanese patients with VSCC to determine the differences in potential somatic mutations among races and identify somatic mutations related to therapeutic targets or prognosis.

Patient characteristics
Nineteen VSCC patients who met the inclusion criteria were enrolled in the NCCH group.The patient characteristics are summarized in Table 1.The age indicates the age at diagnosis.Most of the samples were HPVindependent tumors (84.2%).In this group, all patients with HPV negative tumors had early stage disease.Half of the patients had early-stage disease at the time of diagnosis, whereas the other half had advanced-stage disease.Among the 19 patients, eight experienced recurrence (42.1%) and five died (26.3%).
The clinical characteristics of the 29 patients in the C-CAT group are also summarized in Table 1.The age also indicates the age at diagnosis.Because C-CAT collects data from patients with advanced or progressive disease, most cases were relapsed (96%).
Of the 19 cases, p53 IHC was performed in 16.All three HPV-associated tumors exhibited wild-type staining patterns.In comparison, in 13 cases considered HPV-independent tumors, diffuse strong positive staining was observed in eight cases, and a wild-type pattern was observed in five cases (Supplementary Table 1).
Among 19 cases, three were classified as HPV-associated tumors based on p16 block-type positivity, and 13 were classified as HPV-independent.Regarding preinvasive lesions, classic vulvar intraepithelial neoplasia (VIN) 3 was observed in three cases of HPV-associated tumors.In comparison, differentiated VIN was observed in

Genomic alterations in VSCC detected in the two groups
In the NCCH group, mutation patterns were observed in all HPV-independent tumors, whereas no mutations were found in HPV-dependent tumors (Fig. 1).The gene with the highest frequency of somatic mutations was TP53 (81%), followed by HRAS (26%), CDKN2A (21%), and PIK3CA (5%) (Fig. 1 and Supplementary Table 2).The rate of TP53 mutations in recurrent cases was remarkably high (89%).The candidate somatic mutations in the two groups are listed in Supplementary Tables 2 and 3.In both groups, most TP53 mutations occurred at the DNA-binding domain (DBD), a 191-amino-acid protease-resistant fragment (residues 102-292) containing hotspot mutations at R175, R248, R273, and R282.Most TP53 mutations are located in the DBD 25 , and the mutation profile was similar to those previously reported for VSCC.
In the C-CAT group, TP53 was also the most frequently mutated gene (52%), followed by CDKN2A (24%), PIK3CA (10%), and HRAS (7%) (Fig. 2 and Supplementary Table 3).In the C-CAT and NCCH groups, TP53 mutations were detected in > 50% of patients, showing a trend toward a higher mutation proportion than previously reported (Table 2).The frequency and type of mutations in other genes were similar to those observed in Western populations.Regarding copy number alterations, PIK3CA, SOX2, BLC6, and MAP3K13 amplification was detected in three patients (10%) and CCND1 alterations were observed in two patients (7%).CDKN2A loss was detected in two different cases (7%).

TP53 mutation as a prognostic factor
In the NCCH group, the frequency of TP53 mutations was 8/9 in patients with recurrent disease and 5/9 in patients without recurrent disease.TP53 mutation correlated with a worse prognosis regarding relapse-free survival (log-rank test, P = 0.089) (Fig. 3A), and this association was also observed in patients without HPV infection (log-rank test, P = 0.090) (Fig. 3B).However, did not correlate with worse overall survival in the NCCH and C-CAT groups (log-rank test, P = 0.45, P = 0.13) (Fig. S3).

Discussion
In this study, we investigated the genomic mutational profiles of 48 Japanese VSCC patients and compared them with those of Caucasian patients with VSCC data obtained from previous genomic analyses.The mutational pattern was comparable between Japanese and Caucasian VSCC patients.TP53 was the most frequently mutated gene, followed by HRAS, CDKN2A, and PIK3CA.TP53 mutation correlated with clinical outcome, which supports previous studies identifying TP53 mutation as a prognostic factor.
The frequency of common gene mutations in the present groups did not deviate significantly from that reported previously, except for TP53 mutations, which were observed more frequently.TP53, a crucial tumor suppressor gene, is the most frequently mutated gene in human cancers 26 .Corey et al. detected TP53 mutations in approximately 80% of VSCCs without HPV infection 24 , which is consistent with the present findings.Racial disparity regarding the frequency of TP53 mutations in VSCC may be due to differences in the prevalence of HPV infection.Kang et al. 7 suggested that there is geographical variation in HPV prevalence in vulvar carcinoma.The HPV infection rate in the present group was 16%, whereas that reported previously is approximately 40% 8 .Therefore, the higher frequency of TP53 mutations in this group may be due to the low proportion of HPVassociated VSCC, which is characterized by a low frequency of TP53 mutation.The present findings suggest that genetic factors alone do not account for the regional variation in the epidemiology of vulvar carcinoma.
TP53 mutation is significantly associated with advanced disease and a poor prognosis for various cancer types 27,28 .In this group, the recurrence rate was notably higher in cases with TP53 mutations.Although the number of VSCC patients analyzed was not sufficient to establish TP53 alteration as a poor prognostic factor, the correlation between TP53 alterations and clinical outcomes in VSCC was examined in previous studies (Table 3) [18][19][20][21][22] .In previous studies, TP53 alterations were detected used Sanger sequencing, targeted panel sequencing, and IHC.Nooij et al. 18 reported worse recurrence-free survival in an HPV negative/p53 mutated group (P = 0.044) when comparing three groups based on TP53 mutational status and HPV status (i.e., HPV positive, HPV negative/p53 wildtype, and HPV negative/p53 mutated groups).Kortekaas et al. 19 and Lérias et al. 20 also showed significant differences in recurrence-free survival (P < 0.001) and rate of recurrence (P = 0.0005) according to TP53 mutation status.By contrast, Tessier-Cloutier et al. 22 showed that TP53 mutation alone was not an independent prognostic factor (P = 0.23 for OS and = 0.35 for PFS), whereas TP53 and PIK3CA double mutation was associated with a poor prognosis (P = 0.0016 for OS and P = 0.034 for PFS).Choschzick et al. 21id not identify TP53 as an independent prognostic factor (P > 0.05 for cumulative survival) when comparing three groups with different TP53 mutation and HPV status.In brief, although TP53 mutation correlates with worse prognosis in several studies [19][20][21] , other studies report opposite findings 22,23 .TP53 mutation may be a poor prognostic factor for VSCC, although further studies are needed to confirm the role of TP53 mutation as a critical predictor of worse clinical outcomes.In VSCC, TP53 was the most frequently mutated gene, followed by HRAS, CDKN2A, and PIK3CA, regardless of ethnicity.TP53 mutations, which are present in many cancer types, are considered targets for novel cancer therapies [28][29][30] .A recent study by Guiley 31 reported compounds that directly target mutant p53; however, further studies will be required to determine whether these compounds are therapeutically effective.HRAS, a member of the RAS family of genes, is also known as the mouse sarcoma virus and is associated with human cancer.Tipifarnib is a farnesyltransferase inhibitor that suppresses HRAS function and is used as a targeted therapy for HRAS-mutated tumors 26,32,33 .Ho et al. investigated the efficacy of tipifarnib in head and neck squamous cell carcinoma patients with a mutation in HRAS 34 , reporting an objective response rate of 55%.Loss of the tumor suppressor gene CDKN2A, which encodes p16 and p14, is a frequent occurrence in VSCC.Hsu et al. reported that CDK4/6 inhibitors are a breakthrough therapy for a subset of sarcomas 35 .Therapeutic strategies targeting CDKN2A loss hold great potential for the treatment of various cancer types including VSCC 36 .Pharmacological inhibition of the p16 targets CDK4/6 is a prime example of such a strategy 37 , although additional studies are needed prior to its clinical application.BYL-719 is a PIK3CA inhibitor with low drug-drug interactions, suggesting its potential for combinatorial therapy of recurrent/advanced cervical cancer 38 .Alpelisib (BYL-719) plus fulvestrant was recently approved for the treatment of ER-positive breast cancer patients who develop resistance to ER-targeted therapy 39 .An expansion in the indications for molecular-targeted drugs for VSCC is expected.
This study had several limitations.The number of patients was limited because of the rarity of this cancer.Some clinical data in the C-CAT group were not available, and further validation studies are needed.Additionally, because genomic testing requires a certain amount of tissue, early-stage disease cases and small tumors were not included in the analysis, which may have led to selection bias in the NCCH group.
In conclusion, we found similarities in the somatic mutation profiles of VSCC between Japanese and Caucasian populations.This suggests that the genetic etiology of VSCC does not differ significantly according to ethnic  www.nature.com/scientificreports/background.Analysis of the correlation between TP53 mutation status and clinical outcomes confirmed previous findings identifying TP53 mutation as a prognostic factor.Molecular-targeted therapies may be indicated for selected VSCC patients based on somatic mutational profiles.This study, which included a specific number of Asian patients, will contribute to our understanding of this less investigated yet clinically crucial disease.

Characteristics of 48 patients with VSCC (divided into two groups)
The clinical data used in this study were obtained retrospectively from patient records.The study included 48 Japanese patients with primary VSCC who were divided into two independent groups.The decision trees for patient selection are provided in Fig. S1

Immunohistochemistry (IHC) of p16
IHC was performed using p16 as a surrogate marker to identify HPV-associated tumors.H&E (Fig. S2A and  C) and IHC staining were performed on representative whole sections of tumor from formalin-fixed, paraffinembedded (FFPE) specimens, and the paraffin blocks were sectioned at a thickness of 4 μm.IHC was performed using a Ventana Benchmark Ultra automated immunostainer (Tucson, AZ, USA) and a mouse monoclonal antibody specific for p16 (E6H4, prediluted, Roche diagnostics) according to standard protocols and with appropriate positive and negative controls.A case was classified as an HPV-associated tumor only if it showed blockpositive p16 staining, with diffuse nuclear or nuclear and cytoplasmic positivity (Fig. S2B).Negative results (Fig. S2D) were defined as total absence of staining or weak, focal, and discontinuous staining, which indicated HPV-independent tumors.
The results of p53 immunohistochemical staining analyses were available for 16 out of 19 cases.Staining was performed using a p53 monoclonal antibody (DO7, prediluted, Agilent/Dako, Glostrup, Denmark) in all cases.

DNA preparation and targeted sequencing
FFPE tissue blocks from a VSCC sample were retrieved and serially sectioned at 10 μm thickness on slides.The slides were deparaffinized, and the H&E-stained sections were used for macrodissection.Genomic DNA was extracted from 19 VSCC tissue samples using the QIAamp DNA FFPE tissue kit according to the manufacturer's instructions (Qiagen, Hilden, Germany).Genomic DNA (50 ng) purified from tumor tissues was used for library construction using the Ion AmpliSeq™ Cancer Hotspot Panel v2 (Thermo Fisher Scientific, Waltham, MA, USA).Sequencing was performed on the Ion Proton platform (Thermo Fisher Scientific).For quality control, samples with a mean read depth of coverage > 1000 and a base quality score of 20 (with < 1% probability of being incorrect), which accounted for 90% of total reads, were selected [41][42][43][44] .

HPV genotyping by Sanger sequencing
The 19 samples were subjected to HPV genotyping.Genomic DNA (10 ng) was amplified by PCR targeting two distinct HPV genomic regions using TaKaRa Taq DNA polymerase (Takara Bio Inc., Shiga, Japan).The PCR reaction was performed as described previously 41 .Sanger sequencing was performed using an ABI 3130xl DNA Sequencer (Applied Biosystems).The similarity between the obtained sequences and various HPV genotypes in the GenBank database was determined using the Basic Local Alignment Search Tool (BLAST) (https:// blast.ncbi.nlm.nih.gov/ Blast.cgi).

Figure 1 .
Figure 1.Clinicopathological characteristics and somatic mutations in the NCCH group.The 19 patients are categorized according to clinicopathological features and major somatic mutated genes.Mutated genes are color-coded according to mutation type.

Figure 2 .
Figure 2. Clinical features and genomic alterations in the C-CAT group.The 29 patients are categorized according to available clinical data and major somatic mutated genes.Mutated genes are color-coded according to mutation type.

Table 1 .
Characteristics of patients with vulvar squamous cell carcinoma.NCCH, National Cancer Center Hospital; C-CAT, Center for Cancer Genomics and Advanced Therapeutics; HPV, human papilloma virus; IHC, Immunohistochemistry; PCR, polymerase chain reaction; FIGO, The International Federation of Gynecology and Obstetrics; CCRT, concurrent chemoradiation therapy.The age indicates the age at diagnosis.Lymphovascular invasion was observed in 12 of 18 cases.Lymph node dissection was performed in 16 of the 19 cases, with a median number of excised lymph nodes of 14 (range, 1-56).Lymph node metastasis was observed in 10 of 17 evaluable cases, with three cases showing extracapsular invasion of the lymph nodes.Surgical margins were evaluated in 18 primary tumor resection cases, excluding one lymph node metastasis resection, with 13 cases showing negative margins, five positive horizontal margins, and one positive vertical margin (Supplementary Table1).

Table 2 .
Frequency of mutations in top mutated genes from NCCH, C-CAT, and previous studies.NCCH, National Cancer Center Hospital; C-CAT, Center for Cancer Genomics and Advanced Therapeutics; N/A, not applicable.

Table 3 .
Correlation between TP53 status and clinical outcomes in previous vulvar squamous cell carcinoma studies.NGS, next generation sequencing; IHC, immunohistochemistry. Worst prognosis group is shown in bold.
Vol:.(1234567890) Scientific Reports | (2024) 14:13058 | https://doi.org/10.1038/s41598-024-63913-z 40Forty-two vulvar carcinoma patients were treated at the National Cancer Center Hospital (NCCH) between January 2015 and January 2023.The tumors were diagnosed pathologically according to the 2021 International Federation of Gynecology and Obstetrics criteria, and classified according to the World Health Organization classification of tumors.Of the 42 patients, 23 were excluded due to the histological type or insufficient material; therefore, 19 patients were eligible for inclusion in the NCCH group.Among the 91 patients who were diagnosed between 2019 and 2022 and registered at the Center for Cancer Genomics and Advanced Therapeutics (C-CAT) databank (ver.20221218),30patients with vaginal cancer and 32 patients with other histological types were excluded.Ultimately, 29 patients with VSCC were included in the C-CAT group.In this study, there was no overlap between patients in the NCCH group and those in the C-CAT group.C-CAT is a national data center for cancer genomic medicine in Japan that was established in 201940.Comprehensive genomic profiling tests used were FoundationOne ® CDx and OncoGuide™ NCC Oncopanel.As of June 30, 2022, C-CAT has aggregated 36,340 data points of clinical/genomic information, and is open for sharing with academic institutions and industries.This study was performed in accordance with the Declaration of Helsinki and approved by the institutional review board of the National Cancer Center Research Institute (approval numbers 2017-136, 2020-067) and the C-CAT Data Utilization Review Board (approval number CDU2021-001N).All participants provided written informed consent.